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Abstract: We study oscillons in D + l space-time dimensions using a spherically symmet- 
ric ansatz. From Gaussian initial conditions, these evolve by emitting radiation, approach- 
ing "quasi-breathers" , near-periodic solutions to the equations of motion. Using a truncated 
mode expansion, we numerically determine these quasi-breather solutions in 2 < D < 6 
and the energy dependence on the oscillation frequency. In particular, this energy has a 
minimum, which in turn depends on the number of spatial dimensions. We study the time 
evolution and lifetimes of the resulting quasi-breathers, and show how generic oscillons 
decay into these before disappearing altogether. We comment on the apparent absence of 
oscillons for D > 5 and the possibility of stable solutions for D < 2. 
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1. Introduction 

Topological defects such as monopoles, vortex strings and domains walls owe their stability 
to non-trivial boundary conditions or, equivalently, a conserved (topological) charge Q. 
There are other mechanisms which can lead to the formation of lumps in field theory, 
such as the conservation of a Nother charge allowing for Q-balls or it may simply be 
energetically favourable as in the case of semi-local strings || . 

However, it turns out that even in the absence of such topological constraints or sym- 
metry arguments, solutions to non-linear field equations can have localised lumps of energy. 
A typical example of this is the solitonic breather solution in the sine-Gordon model. 

A more curious set of lumps in field theories are those which exist for a long, but finite, 
time without any of the reasons listed above. These go by the name of oscillons Q and exist 
as a localized oscillating state, with the oscillations remaining coherent for much longer 
than naive dimensional analysis would suggest |5|, Simulations of (ft 4 scalar theories in 
two dimensions @ suggest that they decay only after millions of natural time units whereas 
in three dimensions they last for thousands of natural time units ||. The mechanism for 
their demise follows from their oscillating nature, with the emitted radiation eventually 
causing oscillons to decay. There is precedence for oscillating solutions to continue for 
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ever, although it has only been observed in 1+1 dimensions in breathers of the sine-Gordon 
model. However, this model is rather special in that it is integrable, containing an infinite 
number of conserved currents. An understanding of the evolution of oscillons was gained 
in H using Fourier analysis to construct strictly periodic solutions, later termed "quasi- 
breathers" in [10]. Such exactly periodic solutions exist because they include a component 
of inwardly directed radiation. So while quasi-breathers are not expected to be physical 
their construction shows that the amount of radiation involved is small enough (in some 
cases) to make them a valid approximation. 

In |ll], [l^] it has been suggested that oscillons cease to exist, in any meaningful sense, 
for dimensions greater than 5 or 6 depending on the details of the potential. This result 
was found by assuming a Gaussian profile for oscillons and substituting this into the action 
in order to generate an effective action. Given the extreme sensitivity of oscillon lifetimes 



on the initial profile M, 12] this approach can only yield an indication of what may happen 



and serves as impetus for our study. We use quasi-breathers in our analysis to show how 
oscillons behave in various dimensions, including their non-existence in higher dimensions. 

It is worthwhile mentioning at this point that oscillons have been observed in the 
laboratory [jHj, forming in granular materials placed on a vibrating plate. While these 
oscillons are not described by a relativistic field theory they are similar enough to bear the 
same name. 

Apart from the challenge of understanding oscillons in themselves, interest derives 
from the possible impact of these objects if they are generated in large numbers in phase 
transitions. In particular, if such "coherently oscillating" lumps persist for hundreds or 
thousands of oscillations, they may play a role during post-inflationary preheating and 
baryogenesis, and delay the subsequent equilibration and thermalisation. 

Oscillons have been extensively studied in recent years, mostly in scalar models, but 
an oscillon has also been reported in the SU(2)-Higgs model fl5|| . The very recent work of 
|1C] is a detailed numerical study of oscillons in 3+1 dimensions. Improving on the analysis 
of [}13[ |, the authors propose that generic oscillons are composed of a core quasi-breather of 
a particular oscillation frequency uj and a number of additional modes, depending on the 
initial state. As time goes on the oscillon will radiate the excess energy while changing 
frequency ||, its core scanning through a series of quasi-breathers, eventually approaching 
a critical uj quasi-breather, after which the oscillon decays. The critical uj corresponds to 
the lowest energy quasi-breather. 

The study in two dimensions of Q showed that a non-spherical perturbation was not 
deleterious to the oscillon's evolution and so for our study we restrict our attention to 
spherical symmetry in D space dimensions. The spatial dimension then simply appears 
as a parameter in an effective 1+1 spacetime field theory, this allows us to consider non- 
integer dimensions as a way of observing more clearly how oscillon behaviour depends on 
D. Using this ansatz, we will perform a mode-decomposition of the quasi-breather in a 



way similar to |J, 1C, 13], and numerically find solutions for each u, for each choice of 



spatial dimensions D. Restricting to 3+1 dimensions, it was observed in j|] that there is 
a particular value of qusi-breather frequency, w cr it, which has minimum energy. Oscillons 
with frequencies different from uj CT it evolve toward this solution as they radiate energy, 
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decaying once they reach this frequency. We shall see that this behaviour persists in other 
dimensions, although the time it takes for oscillons to reach u; cr i t depends sensitively on 
D. We shall confirm the fact that for D = 3 lifetimes are of the order of 10 4 [[R]], whereas 
for D = 2 oscillons can oscillate for many times 10 6 . Obviously, in the limit of D = an 
oscillon is simply a particle oscillating in a well and so is strictly periodic. 

These quasi-breather profiles can then be used as initial conditions for real-time sim- 
ulations in order to find the lifetime of these objects. Also, we will use (non fine-tuned) 



Gaussian initial conditions as in [10, 13] and study how these oscillons approach the quasi- 
breather evolution as radiation is emitted. We will see that the oscillon closely follows the 
quasi-breather, and that the two decay in a similar way. 

In section |2| we will set up the model, the quasi-breather ansatz and the mode decom- 
position we use. In section || we will determine the quasi-breather solutions numerically, 
for a range of to and D. The real-time evolution of these quasi-breathers is presented in 
section ||, where we also compare to the evolution of a generic Gaussian initial condition. 
We will discuss the lifetimes of these objects in section |5] and conclude in section ||. 



2. Model, ansatz and mode expansion 

We shall study oscillons in scalar </> 4 in D + 1 dimensions starting with the action 



didPx 



(2.1) 



From this staring point we restrict the number of free parameters by making the potential 
a symmetric double well, achieved by imposing 
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where ^4_d_i = 2ir D / 2 T(D /2) is the volume of the D — 1 sphere which has been integrated 
out. 

A quasi-breather is a periodic solution to the equations of motion, and as such can be 
expanded in modes with frequencies that are integer multiples of some basic frequency u>, 

oo 

(r./j —o n (r)+> d) n (r)cos(nut). (2.6) 
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The in the first term is purely for convenience later. By integrating the action over a 
period, 2it/lj we have 
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We now have an effective action in terms of dimensionless parameters r,u for the dimen- 
sionless fields <j) n . This allows us to reinterpret it in terms of a particle rolling in "time" r 
through a space spanned by the coordinates under the influence of the potential (2.8). A 
contour plot of this potential, setting 02 = 03 = ... = and taking \ < uj < 1 is shown in 
Fig. H, with a maximum at (0o,0i) = (1,~ 1.3), minimum at (1,0) and a saddle at (0,0). 
We shall use this plot to undertand the basic shape of the quasi breather profile. 



3. Finding the quasi-breather pro- 
file 

To find the quasi-breather profile for a par- 
ticular choice of frequency, uj, we need to 
solve for the trajectory of a particle mov- 
ing in U e ff. Note however that, because 
of the r D ~ l in S, there is a frictional force 
acting to slow down the particle in analogy 
to finding the bounce profile in quantum 
tunnelling [16]. Asymptotically, r — > 00, 
we require that the particle takes the value 
= 0, corresponding to the scalar field 
living in the vacuum asymptotically. If we 
were to ignore those n with n > 2 the 
problem would be to find the location in 
Fig. H where a particle can start at "time" 1 




0.6 (S 0.8 



Figure 1: A slice through the effective potential. 
Maxima are at (0,0), (1,1.3) and there is a mini- 
mum at (1,0). 



in order to end at 



as r — > 00. This is 



achieved by taking an initial guess at the starting point and observing whether the particle 
undershoots or overshoots = 0, then amending the starting point appropriately. In the 
actual problem one cannot neglect <p n >2 but the process for finding the values </>o(0) and 
4>i (0) is the same. The starting values of (p n >2 is less critical, as can be understood from 
Q2.8D , and correspond to adjusting the amount of radiation contained in the quasi-breather 
solution. Near = 0, i.e. the asymptotic region, we see that the effective potential is 
quadratic with (f> n aquiring a mass-squared of —5(1 — n 2 ui 2 ), so those n with n > 1/u) 
have positive mass-square while those with n < 1/u; have negative mass-square. This means 



that the 



) n>l/w 



are oscillatory in the asymptotic region, corresponding to radiation. We 



shall only concern ourselves with i < uj < 1 so that n >2 correspond to radiation modes. 
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With (|2,6|) being an expansion of spatially oscillating standing waves, this displays the fact 
that there is ingoing radiation balancing the outgoing. 

3.1 Shooting 



0.8 




Figure 2: Left: The profiles of 0o, 01 1 02 1 4>3 an d 04 f° r D — 3 and w — 0.97. Right: The 
dependence of the quasi-breather radiation on 02 (f = 0) and 03 (r = 0) . 

A general analytical solution of the equation of motion does not exist, and we will 
instead proceed to find the quasi-breather profile by numerical means, using a shooting 
algorithm. Although one can in principle shoot with multiple degrees of freedom we will 
truncate the mode expansion, using the first four terms, <j> n , n =0, 1, 2, 3 in our calculations 
of energy and radiation. Once we have a profile with these components we include the 
fifth field, 04, as a consistency check to make sure that it does indeed form only a small 
contribution. 

In the shooting algorithm we pick initial values for n , n < 4 at r = 0, requiring 
<f>' n (r = 0) = for regularity. As mentioned earlier, the values of the radiation modes, 
n =2, 3, simply dictate how much radiation the quasi breather contains. In [l^] the 
solution is taken which minimizes the amplitude of 02 asymptotically, here we aim to find 
the solution which minimizes the total radiation as discussed below. For any particular 
choice of 02(0) and 03(0) we then vary 0o(O) and 0i(O) until we find the solution which 
approaches = asymptotically. This then gives us the quasi breather for frequency u 
and the chosen 02(0), 03 (0). By finding the quasi breathers for a range of 02 (0), 03 (0) we 
are able to select the one which has the minimum amount of radiation. 

An example of a set of profile functions is shown in Fig. ^ (left); one can see that the 
amplitudes of the high-n modes decreases, allowing for a valid truncation approximation. 
As we vary the values of 02 (r = 0) and 03 (r = 0) in the shooting algorithm we find quasi- 
breathers with differing amounts of radiation, with one choice giving a minimum, Fig. || 
(right). We take this solution with minimum radiation to be our quasi-breather. 

3.2 Radiation energy 

As the oscillon is made by balancing incoming and outgoing radiation then in the realistic 
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situation where there is only outgoing radiation the oscillons will necessarily decay. Here 
we calculate the rate of energy loss as estimated by constructing a quasi-breather. 
First we introduce the standard energy momentum tensor for a real scalar field, 

1 



= dfj,(j)d u (f) - g^ v 

(3.1) 

then we define the momentum flux 

and conservation of the energy-momentum tensor tells us that P^ is divergence free. By 
taking P° as the energy density we find 

E = - J d D x y ^ ) V l P i , (3.3) 

and for a spherically symmetric field we find 

E(R) = -A^R^P^r = R), (3.4) 

where we have defined E(R) to be the energy within a radius R. 

Asymptotically the fields are small so we may approximate the field equation by the 
Klein-Gordon equation 

-4> + cf>" + ^— V = <p, (3.5) 
r 

the closed form solution of which can be approximated asymptotically by the form 

0~ cos(kr-nt + V )/r^ D ^/ 2 , k 2 = Q 2 - 1. (3.6) 

where we have allowed for a possible phase, <p. Our solution is a sum of waves of different 
frequencies producing a standing wave between ingoing and outgoing waves, 

<j)(r -> do) -> (D * 1)/2 (f 2 cos(k 2 r + ip 2 ) cos(2wi) + f 3 cos(k 3 r + ip 3 ) cos(3wt) + ...) (3.7) 

h cos(/c 2 r + ip 2 + 2wt) + f 2 cos(k 2 r + cp 2 - 2ut) (3.8) 



2rP-i)/2 

+/ 3 cos(/c 3 r + <^ 3 + 3wt) + / 3 cos(k 3 r + ip 3 - 3ut) + ...), 
= <fcn + '/'out, (3.9) 
where A; 2 = n 2 w 2 — 1. Using 

P r = kut^'out (3-10) 

we can integrate over one period and find that the average rate of energy loss over a period 
is 

< E > = _^zil {2k2 f 2 + 3k 3 f 2 + ...). (3.11) 
For the profiles used in Fig. [2] we find 

E = 7.9 x 1(T 6 + 0.31 x 1(T 6 + 0.035 x 1(T 6 + ... (3.12) 
The leading term is from the (p 2 , with the (j) 3 and 04 contributions being sub-dominant. 
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3.3 Quasi-breather results 
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Figure 3: The rescaled energy plot. 
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We now have everything in place to find a quasi-breather for each given frequency, we 
shall only focus on frequencies in the range 0.95 < oj < 0.99. The choice of lower bound 
is because, as we shall see, oscillons spend most of their time above this frequency. Our 
upper bound is due to numerical inaccuracies in the shooting algorithm giving unreliable 
results above this frequency. The first result to discuss is the energy of the quasi-breathers 
as u varies, then we can understand how oscillons evolve. 
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Figure 5: The location of the energy minimum. 
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Figure 6: The radiation of a quasi breather. 



In Fig. U we present a plot, for various dimensions, of the quasi-breather energy as a 
function of u. In order to get the data for the different dimensions on the same axes we 
have rescaled each curve so that they meet at u = 0.95. The different rescaling parameters, 
E(D,uj = 0.95), are given in Fig. || showing that to a good approximation the energy of 
a quasi breather depends exponentially on dimension. In order to calculate the energy of 
the quasi-breather we note that the asymptotic form of the field causes the energy, E(R), 
to diverge linearly with R. This divergence is simply due to the radiation, so to define the 
energy of the quasi-breather we simply subtract this contribution. 
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Figure 7: The ratio of radiation to energy as a function of dimension. 



From Fig. || we see that for each dimension considered the energy shows a minimum at 
w cr it as u) is varied, with u; cr it being dimension dependent. This minimum was first observed 
in H and later in |1C]. We can see how this minimum varies with dimension in Fig. || 



Unfortunately our numerical methods break down below D ~ 2.5, because the 4>4 mode 
becomes relevant, and above D ~ 4.5, because we could not shoot (p to the vacuum, so we 
are unable to give a complete plot of the energy minima. It is possible that the techniques 
advocated in [ 10 1 would give a more complete picture, it would be particularly interesting 



to see how the curve develops below D = 2.5. 

As explained in § oscillons are expected to evolve by radiating energy, and so move 
on to a quasi breather with lower energy and different u. Once the quasi-breather with 
w cr it is reached there is no quasi-breather with lower energy, so the oscillon radiation causes 
the lump to radiate away completely. We shall see in the next section that this is a good 
description of the dynamics of oscillons, although we remind the reader that quasi-breathers 
are only expected to be an approximation to oscillon evolution owing to the presence of 
ingoing radiation in a quasi-breather but not in an oscillon. 

Along with the calculation of energy for the quasi-breathers we are able to calculate 
the radiation emitted (and received) as a function of lo for each dimension, this is presented 
in Fig. |6|. A measure of how long an oscillon can be expected to last comes from the ratio 
of the rate of energy loss ( 3.1 1] ) to energy and is plotted in Fig. [7|. From this graph we 



see that as we increase the number of spatial dimensions the amount of radiation relative 
to energy inside a quasi-breather increases by many orders of magnitude, so reducing the 
lifetime of an oscillon. In fact, we can derive a more complete version of the history of 



oscillons using quasi-breathers as we shall see in section 4.1 
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4. Real-time evolution 



Through a detailed minimisation process, we have now determined the solutions with the 
lowest radiation for a range of u and D. The radiation modes correspond to in- and out- 
going radiation, and if they are both included, the quasi-breather should in principle be 
stable, possibly up to the order of truncation of the mode expansion. 

In physical situations we cannot expect there to be incoming radiation so quasi- 
breathers are necessarily an approximation to oscillons. Here we test the state of this 
approximation by comparing the real-time evolution of lumps with initial conditions that 
are either Gaussian or quasi-breather. The quasi-breather profile would, of course, lead 
to an exactly periodic evolution so to make the comparison more sensible we damp the 
radiation at some large radius. 

Ultimately, we are interested in lumps of energy created through some finite energy 
process, such as a phase transition, and so a real-life oscillon is not a quasi-breather. 
Following |]io|| ) we will distinguish between quasi-breathers and oscillons, and represent the 
latter by Gaussian initial conditions. 

We will study the time evolution by discretising the equation of motion, in a simple 
way. The equation of motion reads, imposing spherical symmetry explicitly, 

d t Mt) + eA(Mt)) = n r (t), 

d t u r (t) + eA(n r (t)) = -rn r (t) + d r d' r <t> r (t) - J- (1 - (i - dr/r) D ~ l ) d' r Mt) 

dr ' 



with 



(...)(r,t + dt)-(...)(r,t) , (...)(r,t)-(...)(r,t-dt) 
dt{ - )r - dt ' dt{ - )r ~ dt ' (42) 

and similarly for spatial derivatives, and 

eA(...) =6^ (d r d' r ) 2 (...). (4.3) 

We apply this simple leap-frog algorithm with a large lattice (20001 sites), small lattice 
spacing (dr = 0.01), time-step (dt = 0.005) and an equally spaced grid with boundary 



conditions (f)(r = 20001) = 0, d r (j)(r = 0) = 0. In [ 10 1 it was found that at the level of 
fine-tuning applied there, the maximum lifetime and the corresponding initial condition 
depend somewhat on the discretisation. However at the level of precision employed here, 
this will not be important; we are not attempting to fine-tune to the maximum lifetime. 

In order to get rid of the emitted radiation, we mimic an infinitely distant spatial 
boundary by adding a damping term only beyond a certain radius R b . 

T = 0.1, r> R b ; T = 0, r < R b . (4.4) 

The oscillon core roughly stretches to r = 30 and we concluded that R b = 60 and T = 
0.1 were reasonable choices. In order to get rid of spurious radiation production on the 
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lattice scale, we also included a Kreiss-Oliger <9 4 damping term in the equations of motion, 
controlled by a multiplicative parameter e [13, 17|. 

The initial condition for the oscillon is a Gaussian, 



mr = A^\-—y (4.5) 

which is a two-parameter family of profiles, determined by the amplitude A and the width 
a. We will fix A = 2, i.e. the center of the oscillon starts out in the potential minimum at 
(j) = 2, while the asymptotic region is in the minimum at = 0. We will also fix a = 3.3, 
and consider it a generic initial condition. Not all values of a give long-lived oscillons, but 
3.3 is by no means fine-tuned. 

We will compare the evolution of the Gaussian oscillon to the trajectory starting from 
a quasi-breather. 



4.1 Quasi-breathers and oscillons 




t 

Figure 8: The effective oscillation frequency for a D = 3 oscillon with Gaussian initial conditions 
(black, large oscillations), quasi-breather u> — 0.95 initial conditions (red/grey) and the estimate 
Q ( blu e dashed). The orange dotted line represents w C rit- 



During the real-time evolution the oscillon is not stricly periodic, nevertheless we can 
define an effective frequency of the oscillation by studying the core of the oscillon and 
defining 

= T7 , (4-6) 

"^crossing 

where 5t CTOSS i ng is simply the time between one crossing of 4>(r = 0) through zero and the 
second crossing after that. This will enable us to measure how the frequency of an oscillon 
changes in time, as well as finding the frequency at which the oscillon finally decays. Both 
of these quantities can also be calculated within the quasi-breather picture. The time 
dependence of lo is found because given the energy E(lo) (Fig. |||) we can calculate dE/dw, 
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t t 

Figure 9: The envelope of the center (r = 0) oscillation (left) and the energy within a shell of 
r = Rb = 60 (right). Black lines is the oscillon with Gaussian initial conditions, red/grey using the 
lo = 0.95 quasi-brcather as an initial condition. 

so along with ( |3.11| ) we find t(u) to be 

Fig. |8] presents a comparison in three spatial dimensions of u(t) for our three cases: the 
semi-analytic prediction of |4.7| (blue, dashed); the numerical evolution starting with (non 
fine tuned) Gaussian initial conditions (black); and the numerical evolution starting with 
the Ui = 0.95 quasi-breather profile (red/grey). 

The frequency w e // for the quasi-breather initial conditions slowly scans through uj 
until it reaches u; cr it ~ 0.967 where it decays; this happens around t = 6150. Eq. (|4.7| ) 
predicts a longer lifetime by a factor of ~ 3 giving t ~ 17000, this over-estimation turns out 
to be typical and presumably indicates that the quasi-breather picture underestimates the 
amount of radiation. After the decay, the characteristic frequency becomes 1, the radiation 
frequency. 

The oscillon from Gaussian initial conditions starts out far away from the quasi- 
breather (the center is initially at cfi = 2), but proceeds to shed much of its energy and 
approach the quasi-breather behaviour. The frequency has a large modulation, but as 
the critical frequency is reached the oscillon decays in a way remarkably similar to the 
quasi-breather. 

This approach of the case with Gaussian initial conditions to that of the quasi-breather 
evolution can also be seen by measuring the time development of the core in the lumps, 
Fig. H (left) shows the envelope of the time evolution for the core of the quasi-breather 
(red/grey) and the Gaussian oscillon (black). The modulation of the oscillon frequency is 
manifest, but once the decay happens, the two evolve very similarly. 

Finally, we plot in Fig. ^ (right) the energy within a sphere of radius R b , for the 
Gaussian oscillon (black) and the quasi-breather (red/grey). As expected, the oscillon 
sheds its energy out of the ball into radiation (which is then damped away outside the 
ball). The energy approaches the quasi-breather value, and reaches it just in time for the 
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decay. Note that the decay is delayed by about 600 in time, compared to the crossing of 
the critical frequency. 

Recalling that the initial conditions for the Gaussian were not fine-tuned, this is strong 
evidence that the quasi-breather behaves as an attractor and is a useful way of understand- 
ing oscillon dynamics. 

4.2 More than 3 dimensions 
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Figure 10: Effective frequency (left) and center envelope (right), for D = 4. Color coding as in 
Fig. | 



As argued in 11, 12] and indicated in Fig. [7] we should expect that oscillons have 
rather shorter lives in higher dimensions. The impact of this on the quasi-breather picture 
is that we should expect it to be less accurate, after all, the quasi-breathers are strictly 
periodic and have infinite lifetimes. Repeating the analysis of the previous section, but 



now in four spatial dimensions, we find the results presented in Fig. 1C. Again we note 
that the lifetime as predicted by the quasi-breather arguments is longer than the measured 
lifetime, by a factor of ~ 2 (Fig. |l^ (left)). However, by observing the evolution of the core 
(Fig. |l0| (right)) we see that the oscillon phase ends at around t ~ 400 which is where oj e R 
hits w cr it ~ 0.971, thus strengthening the argument that oscillons decay once they reach 
the critical frequency. 



4.3 Less than 3 dimensions 

Because of concerns about the validity of the mode truncation, we did not trust finding 
quasi-breathers in D = 1, 2; the techniques of |l(]] may well be more reliable in this regime. 
Instead we generalised to non-integer D, D = 2.5, 2.75 which allows us to see more clearly 
any dependence on spatial dimension. We are able to do this becuase we have assumed 
spherical symmetry, meaning that the dimension D appears only as a parameter in our 
evolution equations. 



In Fig. |ll] we see the evolution of u> for oscillons with quasi-breather initial condi- 
tions, Gaussian initial conditions and the semi-analytic estimate. As expected we see that 
oscillons last longer in these smaller spatial dimensions, with the semi-analytic reasoning 
overestimating the lifetime by a factor of a few. 
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Figure 11: Effective frequency for D = 2.75 (left) and D = 2.5 (right). Quasi-breather initial con- 
ditions (red/grey), Gaussian initial conditions (black), estimate (blue dashed) and critical frequency 
(orange dotted). 



Fig. 11 (right) shows u e g(t) for D = 2.5. Note that the timescale of interest is now 



10 b and we did not evolve the simulation up to the decay time in this case; w cr ;t — 0.98 is 
far from reached. The Gaussian oscillon (black) closely follows the quasi-breather oscillon 
from times ~ 20000. The estimate of the evolution (blue dashed) is somewhat off in the 
same way as the other cases, although the timescale and general trend is well reproduced. 
The prediction is that the oscillons could live as long as 8 x 10 6 in natural time units. 

4.4 Approaching u; cr it from above 
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Figure 12: The effective frequency for various D when starting from an uj = 0.99 quasi-breather. 
The orange dotted line is w C rit- 



So far we have considered only those oscillons whose frequency is less than the critical 
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frequency, and we have seen that as they evolve they follow along one of the curves in 
Fig. H from some ui < u> cr it until they reach lo = o; cr it at which point they decay. The 
evolution of the oscillon as it approaches w cr it is observed to be rather gentle, slowing down 
because the radiation of a quasi-breather decreases as we approach w cr it from below. Fig. 
|3| reveals another possibility for oscillons, namely those starting from to > cj cr it, in this case 
the quasi-breather picture would say that u approaches w cr it from above with the oscillon 
following down one of the curves in Fig. ||. Using Gaussian initial conditions we were 
not able to find an oscillon with frequency larger than the critical value, however we are 
able to start with a quasi-breather profile satisfying u > u> cr it- The results are somewhat 
disappointing for the quasi-breather picture, as the subsequent evolution did not show any 
oscillon phase. Instead, with the exception of D = 3, the field dissipated rapidly straight 
into radiation, to = 1, as shown in Fig. |l^. 



5. Oscillon life-times and dimensionality 
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Figure 13: Left: The evolution of u> as predicted by the quasi-breather model using 4.7, Right: 
Lifetimes for an oscillon starting from a Gaussian with a = 3.3 for 2.5 < D < 7. 

In this section we shall collect together the information we have found about oscillons 
and quasi-breathers, and show how the quasi-breather picture predicts oscillon evolution 
should depend on dimension. 

We have already seen some specific examples of how the quasi-breather model predicts 
the evolution of uj using ( |4.7| ), in Fig. [l3| (left) we collect together the predictions for a 
range of dimensions 2^ < D > 5. As noted previously these tend to be off by a factor 
of a few, overestimating the lifetime, but we can clearly see that the oscillons in lower 
dimensions take much longer to reach the critical frequency. In higher dimensions we find 
that the predicted lifetime is simply not long enough for the system to enter into an oscillon 
phase. 

We may now compare the lifetime for some specific initial condition to that of the 
quasi-breather oscillon model and see how they compare in different dimensions 



As was shown in [10, in D=3 one can tune a to produce oscillons with a large 
lifetime, but the level of fine-tuning necessary is extreme. It is then surprising to find that 
in D < 3, oscillons living for t ~ 10 6 are rather more generic and require little tuning of 
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initial conditions. The converse is true in higher dimensions, with the amount of precision 
required to find a long-lived oscillon increasing with D. This is the main reason that the 
results of 12|, based on substituting an oscillon profile into the action, needed to be 



confirmed with the more detailed analysis we have presented. Rather than trying to find 
the fine-tuned initial conditions which maximize the lifetime in each dimension, we shall 
take the same Gaussian initial conditions {a = 3.3, A = 2) in the various dimensions and 
compare the lifetime of the subsequent oscillon with that of the quasi-breather model. The 
results are plotted in Fig. [l3| (right) and show that the semi-analytic predictions are well 
matched by the numerical evolution. 

In order to guide the eye we have overlaid a curve of the form 

c 

^lifetime = 77^ TW — a ' V^-l) 

(i^crit - L>) 

with c = 1000, a = — 1.5 and -D cr it = 2.65, suggesting that the lifetime goes to infinity 
for some D > 0, and possibly for D > 2. While this is a possibility, the quasi-breather 
profile could conspire to have zero amplitude for the radiation modes, we believe that the 
extrapolation of the quasi-breather results to D = 2 require a finite lifetime. In particular, 
Fig. [?] does not suggest that the radiation goes to zero at D = 2. Note also that there is 
a kink in the curve of Fig. ^ (right) at D ~ 5, after which the lifetime suddenly drops 
(there are now only a handfull of periods), where a similar fit gives c = 45, a = —0.7 and 
-Dcrit = 5.02. Although the precise fitting parameters are not essential, the change in power 
suggests a different regime (no oscillons) has set in. 



6. Conclusion 



The interpretation of oscillons in terms of quasi-breathers gives a simple and consistent 
picture of why such long-lived objects should exist and may even be generated in a generic 
phase transition. 

In our setup, the quasi-breathers are the periodic solutions to the equations of motion 
with minimal outgoing radiation. The presence of this radiation renders the quasi-breather 
unstable, although since the radiation is minimal, presumably it is the longest lived. Even 
without excessive fine-tuning, localised lumps (oscillons) can approach and follow the evo- 
lution of quasi-breathers, which hence must be considered (at least) weak attractors in the 
space of classical solutions. 

A quasi-breather exists for each choice of frequency. Given an initial choice of lo, 
the configuration evolves through a sequence of quasi-breathers, approaching a critical 
frequency, corresponding to the lowest energy quasi-breather. The time is takes to do 
so is well reproduced by considering the emitted radiation. As the critical frequency is 
reached, the energy can no longer be lowered, and the quasi-breather decays completely 
into radiation. 

We performed this analysis for a range of D, the number of spatial dimensions, in 
order to establish the stability/existence of oscillons in higher dimensions. Using the quasi- 
breather approach we have provided an understanding for the prescence of very long-lived 
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solutions in low dimensions and short-lived solutions in higher dimensions. Although our 
method could not cope with D < 2, it is tempting to conjecture from (|5.1| ) that the lifetime 
will diverge for some D > 0, although the extrapolation of Fig. [?] does not seem to support 
this. Perhaps including more radiation modes and more hard work can settle this. While 
we have studied the dependence on D we have not studied the dependence on V{4>). One 
especially interesting case would be the sine-Gordon potential where we know there exist 
exactly periodic solutions in D = 1 [pl| as well as D = 0, which is just a particle in a well. 

For D > 3, the lifetimes become smaller and smaller, and oscillons rush to reach the 
quasi-breather before decaying. For D > 5 we were unable to find quasi-breathers with 
our shooting method, and it indeed seems like a wholly new, oscillon-free, regime sets in. 
This is where the time a configuration takes to reach a quasi-breather is longer than the 
lifetime of the quasi-breather itself, so no oscillon can get established. 

In a 3+1-dimensional expanding Universe, assuming that the expansion does not sig- 
nificantly influence the field evolution, it is therefore natural to expect oscillons to be 
generated in phase transitions or during (p)reheating after inflation, with lifetimes of order 
10 in mass units, but probably not much longer. Scalar fields thermalise on much longer 



timescales than this [19, 20 1, but when including gauge degrees of freedom equilibraton is 



much faster (see for instance [21, [2^]). Oscillons may potentially delay this process if they 
are produced in large numbers, and if the effects of the gauge field do not substantially 
change the picture. 

In order to make further statements about this point, we need to look for quasi- 
breathers in more realistic settings such as the Standard Model and its extensions. In the 



restriction to the 5C/(2)-Higgs model [15], an approach similar to the one employed here 
should work, in spite of the complexity introduced by the additional degrees of freedom. 
Work is under way to address this. 
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1. Introduction 

Topological defects such as monopoles, vortex strings and domains walls owe their stability 
to non-trivial boundary conditions or, equivalently, a conserved (topological) charge [?]. 
There are other mechanisms which can lead to the formation of lumps in field theory, 
such as the conservation of a Nother charge allowing for Q-balls [?] or it may simply be 
energetically favourable as in the case of semi-local strings [?]. 

However, it turns out that even in the absence of such topological constraints or sym- 
metry arguments, solutions to non-linear field equations can have localised lumps of energy. 
A typical example of this is the solitonic breather solution in the sine-Gordon model. 

A more curious set of lumps in field theories are those which exist for a long, but finite, 
time without any of the reasons listed above. These go by the name of oscillons [?] and exist 
as a localized oscillating state, with the oscillations remaining coherent for much longer 
than naive dimensional analysis would suggest [?, ?]. Simulations of cf) 4 scalar theories in 
two dimensions [?] suggest that they decay only after millions of natural time units whereas 
in three dimensions they last for thousands of natural time units [?]. The mechanism for 
their demise follows from their oscillating nature, with the emitted radiation eventually 
causing oscillons to decay. There is precedence for oscillating solutions to continue for 
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ever, although it has only been observed in 1+1 dimensions in breathers of the sine-Gordon 
model. However, this model is rather special in that it is integrable, containing an infinite 
number of conserved currents. An understanding of the evolution of oscillons was gained 
in [?] using Fourier analysis to construct strictly periodic solutions, later termed "quasi- 
breathers" in [?]. Such exactly periodic solutions exist because they include a component 
of inwardly directed radiation. So while quasi-breathers are not expected to be physical 
their construction shows that the amount of radiation involved is small enough (in some 
cases) to make them a valid approximation. 

In [?, ?] it has been suggested that oscillons cease to exist, in any meaningful sense, 
for dimensions greater than 5 or 6 depending on the details of the potential. This result 
was found by assuming a Gaussian profile for oscillons and substituting this into the action 
in order to generate an effective action. Given the extreme sensitivity of oscillon lifetimes 
on the initial profile [?, ?] this approach can only yield an indication of what may happen 
and serves as impetus for our study. We use quasi-breathers in our analysis to show how 
oscillons behave in various dimensions, including their non-existence in higher dimensions. 

It is worthwhile mentioning at this point that oscillons have been observed in the 
laboratory [?], forming in granular materials placed on a vibrating plate. While these 
oscillons are not described by a relativistic field theory they are similar enough to bear the 
same name. 

Apart from the challenge of understanding oscillons in themselves, interest derives 
from the possible impact of these objects if they are generated in large numbers in phase 
transitions. In particular, if such "coherently oscillating" lumps persist for hundreds or 
thousands of oscillations, they may play a role during post-inflationary preheating and 
baryogenesis, and delay the subsequent equilibration and thermalisation. 

Oscillons have been extensively studied in recent years, mostly in scalar models, but an 
oscillon has also been reported in the SU(2)-Higgs model [?]. The very recent work of [?] is 
a detailed numerical study of oscillons in 3+1 dimensions. Improving on the analysis of [?], 
they propose that generic oscillons are composed of a core quasi-breather of a particular 
oscillation frequency to and a number of additional modes, depending on the initial state. 
As time goes on the oscillon will radiate the excess energy while changing frequency [?], 
its core scanning through a series of quasi-breathers, eventually approaching a critical uj 
quasi-breather, after which the oscillon decays. The critical uj corresponds to the lowest 
energy quasi-breather. 

The study in two dimensions of [?] showed that a non-spherical perturbation was 
not deleterious to oscillon's evolution and so for our study we restrict our attention to 
spherical symmetry in D space dimensions. The spatial dimension then simply appears 
as a parameter in an effective 1+1 spacetime field theory, this allows us to consider non- 
integer dimensions as a way of observing more clearly how oscillon behaviour depends 
on D. Using this ansatz, we will perform a mode-decomposition of the quasi-breather in 
a way similar to [?, ?, ?], and numerically find solutions for each to, for each choice of 
spatial dimensions D. Restricting to 3+1 dimensions, it was observed in [?] that there is 
a particular value of qusi-breather frequency, uJcrit-, which has minimum energy. Oscillons 
with frequencies different from uj cr it evolve toward this solution as they radiate energy, 
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decaying once they reach this frequency. We shall see that this behaviour persists in other 
dimensions, although the time it takes for oscillons to reach oo C rit depends sensitively on 
D. We shall confirm the fact that for D = 3 lifetimes are of the order of 10 4 [?], whereas 
for D = 2 oscillons can oscillate for many times 10 6 . Obviously, in the limit of D = an 
oscillon is simply a particle oscillating in a well and so is strictly periodic. 

These quasi-breather profiles can then be used as initial conditions for real-time sim- 
ulations in order to find the lifetime of these objects. Also, we will use (non fine-tuned) 
Gaussian initial conditions as in [?, ?] and study how these oscillons approach the quasi- 
breather evolution as radiation is emitted. We will see that the oscillon closely follows the 
quasi-breather, and that the two decay in a similar way. 

In section |2| we will set up the model, the quasi-breather ansatz and the mode decom- 
position we use. In section || we will determine the quasi-breather solutions numerically, 
for a range of u and D. The real-time evolution of these quasi-breathers is presented in 
section ||, where we also compare to the evolution of a generic Gaussian initial condition. 
We will discuss the lifetimes of these objects in section |5| and conclude in section ||. 

2. Model, ansatz and mode expansion 

We shall study oscillons in scalar d) A in D + 1 dimensions starting with the action 

(2.1) 



did D x 



From this staring point we restrict the number of free parameters by making the potential 
a symmetric double well, achieved by imposing 



a = -3mi//3/2. (2.2) 
This gives a potential with extrema at 

<Pext = 0, —=,mJ-. (2.3) 



W V P 

It is now possible to perform some rescalings of the co-ordinates and the field 

A 777, 

d> = — = — <b, r = mf, t = mt, (2.4) 
and consider only radially symmetric field configurations to find 



S = -"HlJLad^ I dt dr r D - 1 



(2.5) 



where ^4_d-i = 2tt d / 2 T(D /2) is the volume of the D — 1 sphere which has been integrated 
out. 

A quasi-breather is a periodic solution to the equations of motion, and as such can be 
expanded in modes with frequencies that are integer multiples of some basic frequency u>, 

1 x 

0(r, t) = -j=Mr) + Yl cos (^)- ( 2 - 6 ) 

n=l 
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The in the first term is purely for convenience later. By integrating the action over a 
period, 2it/lj we have 
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We now have an effective action in terms of dimensionless parameters r,u for the dimen- 
sionless fields <j) n . This allows us to reinterpret it in terms of a particle rolling in "time" r 
through a space spanned by the coordinates under the influence of the potential (2.8). A 
contour plot of this potential, setting 02 = 03 = ... = and taking \ < uj < 1 is shown in 
Fig. H, with a maximum at (0o,0i) = (1,~ 1.3), minimum at (1,0) and a saddle at (0,0). 
We shall use this plot to undertand the basic shape of the quasi breather profile. 




3. Finding the quasi-breather pro- 
file 

To find the quasi-breather profile for a par- 
ticular choice of frequency, uj, we need to 
solve for the trajectory of a particle mov- 
ing in U e ff. Note however that, because 
of the r D ~ l in S, there is a frictional force 
acting to slow down the particle in analogy 
to finding the bounce profile in quantum 
tunnelling [?]. Asymptotically, r — > 00, we 
require that the particle takes the value 
= 0) corresponding to the scalar field 
living in the vacuum asymptotically. If we 
were to ignore those 4> n with n > 2 the 
problem would be to find the location in 
Fig. H where a particle can start at "time" r = in order to end at cj) = as r — > 00. This is 
achieved by taking an initial guess at the starting point and observing whether the particle 
undershoots or overshoots = 0, then amending the starting point appropriately. In the 
actual problem one cannot neglect 4> n >2 but the process for finding the values 4>o(0) and 
4>i (0) is the same. The starting values of (p n >2 is less critical, as can be understood from 
Q2.8D , and correspond to adjusting the amount of radiation contained in the quasi-breather 
solution. Near = 0, i.e. the asymptotic region, we see that the effective potential is 
quadratic with (f> n aquiring a mass-squared of —5(1 — n 2 ui 2 ), so those <fi n with n > 1/u) 
have positive mass-square while those with n < 1/u; have negative mass-square. This means 



0.6 (S 0.8 



Figure 1: A slice through the effective potential. 
Maxima are at (0,0), (1,1.3) and there is a mini- 
mum at (1,0). 



that the 



) n>l/w 



are oscillatory in the asymptotic region, corresponding to radiation. We 



shall only concern ourselves with i < uj < 1 so that 4> n >2 correspond to radiation modes. 
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With (|2,6|) being an expansion of spatially oscillating standing waves, this displays the fact 
that there is ingoing radiation balancing the outgoing. 

3.1 Shooting 



0.8 




Figure 2: Left: The profiles of 4>o, <Pi, 02 1 4>3 an d 04 f° r D — 3 and w — 0.97. Right: The 
dependence of the quasi-breather radiation on 02 (f = 0) and 03 (r = 0) . 

A general analytical solution of the equation of motion does not exist, and we will 
instead proceed to find the quasi-breather profile by numerical means, using a shooting 
algorithm. Although one can in principle shoot with multiple degrees of freedom we will 
truncate the mode expansion, using the first four terms, <j> n , n =0, 1, 2, 3 in our calculations 
of energy and radiation. Once we have a profile with these components we include the 
fifth field, 04, as a consistency check to make sure that it does indeed form only a small 
contribution. 

In the shooting algorithm we pick initial values for n , n < 4 at r = 0, requiring 
<j)' n (r = 0) = for regularity. As mentioned earlier, the values of the radiation modes, 
n =2, 3, simply dictate how much radiation the quasi breather contains. In [?] the solution 
is taken which minimizes the amplitude of 02 asymptotically, here we aim to find the 
solution which minimizes the total radiation as discussed below. For any particular choice 
of 02(0) and 03(0) we then vary 0o(O) and 0i(O) until we find the solution which approaches 
= asymptotically. This then gives us the quasi breather for frequency 00 and the chosen 
02(0), 03(0). By finding the quasi breathers for a range of 02(0), 03(0) we are able to select 
the one which has the minimum amount of radiation. 

An example of a set of profile functions is shown in Fig. ^ (left); one can see that the 
amplitudes of the high-ra modes decreases, allowing for a valid truncation approximation. 
As we vary the values of 02 (r = 0) and 03 (r = 0) in the shooting algorithm we find quasi- 
breathers with differing amounts of radiation, with one choice giving a minimum, Fig. || 
(right). We take this solution with minimum radiation to be our quasi-breather. 

3.2 Radiation energy 

As the oscillon is made by balancing incoming and outgoing radiation then in the realistic 
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situation where there is only outgoing radiation the oscillons will necessarily decay. Here 
we calculate the rate of energy loss as estimated by constructing a quasi-breather. 
First we introduce the standard energy momentum tensor for a real scalar field, 

1 



= dfj,(j)d u (f) - g^ v 

(3.1) 

then we define the momentum flux 

and conservation of the energy-momentum tensor tells us that P^ is divergence free. By 
taking P° as the energy density we find 

E = - J d D x y ^ ) V l P i , (3.3) 

and for a spherically symmetric field we find 

E(R) = -A^R^P^r = R), (3.4) 

where we have defined E(R) to be the energy within a radius R. 

Asymptotically the fields are small so we may approximate the field equation by the 
Klein-Gordon equation 

-4> + cf>" + ^— V = <p, (3.5) 
r 

the closed form solution of which can be approximated asymptotically by the form 

0~ cos(kr-nt + V )/r^ D ^/ 2 , k 2 = Q 2 - 1. (3.6) 

where we have allowed for a possible phase, <p. Our solution is a sum of waves of different 
frequencies producing a standing wave between ingoing and outgoing waves, 

<j)(r -> do) -> (D * 1)/2 (f 2 cos(k 2 r + ip 2 ) cos(2wi) + f 3 cos(k 3 r + ip 3 ) cos(3wt) + ...) (3.7) 

h cos(/c 2 r + ip 2 + 2wt) + f 2 cos(k 2 r + cp 2 - 2ut) (3.8) 



2rP-i)/2 

+/ 3 cos(/c 3 r + <^ 3 + 3wt) + / 3 cos(k 3 r + ip 3 - 3ut) + ...), 
= <fcn + '/'out, (3.9) 
where A; 2 = n 2 w 2 — 1. Using 

P r = kut^'out (3-10) 

we can integrate over one period and find that the average rate of energy loss over a period 
is 

< E > = _^zil {2k2 f 2 + 3k 3 f 2 + ...). (3.11) 
For the profiles used in Fig. [2] we find 

E = 7.9 x 1(T 6 + 0.31 x 1(T 6 + 0.035 x 1(T 6 + ... (3.12) 
The leading term is from the (p 2 , with the (j) 3 and 04 contributions being sub-dominant. 
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3.3 Quasi-breather results 
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Figure 3: The rescaled energy plot. 
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We now have everything in place to find a quasi-breather for each given frequency, we 
shall only focus on frequencies in the range 0.95 < oj < 0.99. The choice of lower bound 
is because, as we shall see, oscillons spend most of their time above this frequency. Our 
upper bound is due to numerical inaccuracies in the shooting algorithm giving unreliable 
results above this frequency. The first result to discuss is the energy of the quasi-breathers 
as u varies, then we can understand how oscillons evolve. 
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Figure 6: The radiation of a quasi breather. 



In Fig. U we present a plot, for various dimensions, of the quasi-breather energy as a 
function of u. In order to get the data for the different dimensions on the same axes we 
have rescaled each curve so that they meet at u = 0.95. The different rescaling parameters, 
E(D,uj = 0.95), are given in Fig. || showing that to a good approximation the energy of 
a quasi breather depends exponentially on dimension. In order to calculate the energy of 
the quasi-breather we note that the asymptotic form of the field causes the energy, E(R), 
to diverge linearly with R. This divergence is simply due to the radiation, so to define the 
energy of the quasi-breather we simply subtract this contribution. 
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Figure 7: The ratio of radiation to energy as a function of dimension. 



From Fig. || we see that for each dimension considered the energy shows a minimum 
at oocrit as u is varied, with uj cr it being dimension dependent. This minimum was first 
observed in [?] and later in [?]. We can see how this minimum varies with dimension in 
Fig. H Unfortunately our numerical methods break down below D ~ 2.5, because the c/>4 
mode becomes relevant, and above D ~ 4.5, because we could not shoot <j) to the vacuum, 
so we are unable to give a complete plot of the energy minima. It is possible that the 
techniques advocated in [?] would give a more complete picture, it would be particularly 
interesting to see how the curve develops below D = 2.5. 

As explained in [?] oscillons are expected to evolve by radiating energy, and so move 
on to a quasi breather with lower energy and different u. Once the quasi-breather with 
Ucrit is reached there is no quasi-breather with lower energy, so the oscillon radiation causes 
the lump to radiate away completely. We shall see in the next section that this is a good 
description of the dynamics of oscillons, although we remind the reader that quasi-breathers 
are only expected to be an approximation to oscillon evolution owing to the presence of 
ingoing radiation in a quasi-breather but not in an oscillon. 

Along with the calculation of energy for the quasi-breathers we are able to calculate 
the radiation emitted (and received) as a function of lo for each dimension, this is presented 
in Fig. |6| A measure of how long an oscillon can be expected to last comes from the ratio 
of the rate of energy loss ( 3.1 1] ) to energy and is plotted in Fig. [7|. From this graph we 



see that as we increase the number of spatial dimensions the amount of radiation relative 
to energy inside a quasi-breather increases by many orders of magnitude, so reducing the 
lifetime of an oscillon. In fact, we can derive a more complete version of the history of 



oscillons using quasi-breathers as we shall see in section 4.1 
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4. Real-time evolution 



Through a detailed minimisation process, we have now determined the solutions with the 
lowest radiation for a range of to and D. The radiation modes correspond to in- and out- 
going radiation, and if they are both included, the quasi-breather should in principle be 
stable, possibly up to the order of truncation of the mode expansion. 

In physical situations we cannot expect there to be incoming radiation so quasi- 
breathers are necessarily an approximation to oscillons. Here we test the state of this 
approximation by comparing the real-time evolution of lumps with initial conditions that 
are either Gaussian or quasi-breather. The quasi-breather profile would, of course, lead 
to an exactly periodic evolution so to make the comparison more sensible we damp the 
radiation at some large radius. 

Ultimately, we are interested in lumps of energy created through some finite energy 
process, such as a phase transition, and so a real-life oscillon is not a quasi-breather. 
Following [?], we will distinguish between quasi-breathers and oscillons, and represent the 
latter by Gaussian initial conditions. 

We will study the time evolution by discretising the equation of motion, in a simple 
way. The equation of motion reads, imposing spherical symmetry explicitly, 

d t Mt) + ^(Mt)) = n r (t), 

d t U r (t) + eA(n r (t)) = -iTI r (t) + d r d' r <f> r (t) ~ ^ (I - (I - dr/r) D ~ l ) d' r 4> r {t) 

- (i - \ut) + 1$®) Mt), (4-1) 



with 

(...)(r,t + dt)-(...)(r,t) , (...)(r,t)-(...)(r,t-dt) 
9t{ - )r - dt ' &t{ - )r - dt ' (42) 

and similarly for spatial derivatives, and 

eA(...) =e^ (d r d' r ) 2 (•••)• (4-3) 

We apply this simple leap-frog algorithm with a large lattice (20001 sites), small lattice 
spacing (dr = 0.01), time-step (dt = 0.005) and an equally spaced grid with boundary 
conditions <fi(r = 20001) = 0, d r (j)(r = 0) = 0. In [?] it was found that at the level of 
fine-tuning applied there, the maximum lifetime and the corresponding initial condition 
depend somewhat on the discretisation. However at the level of precision employed here, 
this will not be important; we are not attempting to fine-tune to the maximum lifetime. 

In order to get rid of the emitted radiation, we mimic an infinitely distant spatial 
boundary by adding a damping term only beyond a certain radius R b . 

T = 0.1, r > R b ; T = 0, r < R b . (4.4) 

The oscillon core roughly stretches to r = 30 and we concluded that R b = 60 and T = 
0.1 were reasonable choices. In order to get rid of spurious radiation production on the 
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lattice scale, we also included a Kreiss-Oliger damping term in the equations of motion, 
controlled by a multiplicative parameter e [?, ?]. 

The initial condition for the oscillon is a Gaussian, 

0(O) r =Aexp(-^Y (4.5) 

which is a two-parameter family of profiles, determined by the amplitude A and the width 
a. We will fix A = 2, i.e. the center of the oscillon starts out in the potential minimum at 
(j) = 2, while the asymptotic region is in the minimum at = 0. We will also fix a = 3.3, 
and consider it a generic initial condition. Not all values of a give long-lived oscillons, but 
3.3 is by no means fine-tuned. 

We will compare the evolution of the Gaussian oscillon to the trajectory starting from 
a quasi-breather. 



4.1 Quasi-breathers and oscillons 




During the real-time evolution the oscillon is not stricly periodic, nevertheless we can 
define an effective frequency of the oscillation by studying the core of the oscillon and 
defining 

Waff(t) = TT ^ , (4-6) 
"'-crossing 

where (^crossing is simply the time between one crossing of <f>{r = 0) through zero and the 
second crossing after that. This will enable us to measure how the frequency of an oscillon 
changes in time, as well as finding the frequency at which the oscillon finally decays. Both 
of these quantities can also be calculated within the quasi-breather picture. The time 
dependence of lo is found because given the energy E(u) (Fig. |||) we can calculate dE/dw, 



-11- 




2000 4000 6000 8000 10000 u 2000 4000 6000 8000 10000 



t t 

Figure 9: The envelope of the center (r = 0) oscillation (left) and the energy within a shell of 
r = Rb = 60 (right). Black lines are the oscillon with Gaussian initial conditions, red the oscillon 
with the uj = 0.95 quasi-breather as an initial condition. 

so along with ( |3.11| ) we find t(uj) to be 

Fig. |8] presents a comparison in three spatial dimensions of uj(t) for our three cases: the 
semi-analytic prediction of |4.7| (blue, dashed); the numerical evolution starting with (non 
fine tuned) Gaussian initial conditions (black); and the numerical evolution starting with 
the LOi = 0.95 quasi-breather profile (red). 

The frequency oj e ff for the quasi-breather initial conditions slowly scans through uj 
until it reaches u; cr it ~ 0.967 where it decays; this happens around t = 6150. Eq. ( |I~^ ) 
predicts a longer lifetime by a factor of ~ 3 giving t ~ 17000, this over-estimation turns out 
to be typical and presumably indicates that the quasi-breather picture underestimates the 
amount of radiation. After the decay, the characteristic frequency becomes 1, the radiation 
frequency. 

The oscillon from Gaussian initial conditions starts out far away from the quasi- 
breather (the center is initially at eft = 2), but proceeds to shed much of its energy and 
approach the quasi-breather behaviour. The frequency has a large modulation, but as 
the critical frequency is reached the oscillon decays in a way remarkably similar to the 
quasi-breather. 

This approach of the case with Gaussian initial conditions to that of the quasi-breather 
evolution can also be seen by measuring the time development of the core in the lumps, Fig. 
|9] (left) shows the envelope of the time evolution for the core of the quasi-breather (red) 
and the Gaussian oscillon (black). The modulation of the oscillon frequency is manifest, 
but once the decay happens, the two evolve very similarly. 

Finally, we plot in Fig. ^ (right) the energy within a sphere of radius Rf,, for the 
Gaussian oscillon (black) and the quasi-breather oscillon (red). As expected, the oscillon 
sheds its energy out of the ball into radiation (which is then damped away outside the 
ball). The energy approaches the quasi-breather value, and reaches it just in time for the 
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decay. Note that the decay is delayed by about 600 in time, compared to the crossing of 
the critical frequency. 

Recalling that the initial conditions for the Gaussian were not fine-tuned, this is strong 
evidence that the quasi-breather behaves as an attractor and is a useful way of understand- 
ing oscillon dynamics. 

4.2 More than 3 dimensions 
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Figure 10: Effective frequency (left) and center envelope (right), for D = 4. Color coding as in 
Fig. | 

As argued in [?, ?] and indicated in Fig. ^ we should expect that oscillons have rather 
shorter lives in higher dimensions. The impact of this on the quasi-breather picture is that 
we should expect it to be less accurate, after all, the quasi-breathers are strictly periodic 
and have infinite lifetimes. Repeating the analysis of the previous section, but now in four 



spatial dimensions, we find the results presented in Fig. 1C. Again we note that the lifetime 
as predicted by the quasi-breather arguments is longer than the measured lifetime, by a 
factor of ~ 2 (Fig. |l0| (left)). However, by observing the evolution of the core (Fig. |l~0| 
(right)) we see that the oscillon phase ends at around t ~ 400 which is where u e ff bits 
^crit ~ 0.971, thus strengthening the argument that oscillons decay once they reach the 
critical frequency. 



4.3 Less than 3 dimensions 

Because of concerns about the validity of the mode truncation, we did not trust finding 
quasi-breathers in D = 1, 2; the techniques of [?] may well be more reliable in this regime. 
Instead we generalised to non-integer D, D = 2.5, 2.75 which allows us to see more clearly 
any dependence on spatial dimension. We are able to do this becuase we have assumed 
spherical symmetry, meaning that the dimension D appears only as a parameter in our 
evolution equations. 



In Fig. |ll] we see the evolution of u> for oscillons with quasi-breather initial condi- 
tions, Gaussian initial conditions and the semi-analytic estimate. As expected we see that 
oscillons last longer in these smaller spatial dimensions, with the semi-analytic reasoning 
overestimating the lifetime by a factor of a few. 
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Figure 11: Effective frequency for D = 2.75 (left) and D = 2.5 (right). Quasi-breather initial 
conditions (red), Gaussian initial conditions (black), estimate (blue) and critical frequency (orange). 



Fig. |IT] (right) shows uj e g(t) for D = 2.5. Note that the timescale of interest is now 
10 6 and we did not evolve the simulation up to the decay time in this case; w cr ;t — 0.98 is 
far from reached. The Gaussian oscillon (black) closely follows the quasi-breather oscillon 
from times ~ 20000. the estimate of the evolution (blue) is somewhat off in the same 
way as the other cases, although the timescale and general trend is well reproduced. The 
prediction is that the oscillons could live as long as 8 x 10 6 in natural time units. 

4.4 Approaching u; cr it from above 
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Figure 12: The effective frequency for various D when starting from an ui = 0.99 quasi-breather. 
The orange dotted line is w C rit- 



So far we have considered only those oscillons whose frequency is less than the critical 
frequency, and we have seen that as they evolve they follow along one of the curves in 
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Fig. [3| from some u < uj cr it until they reach u> = u> C rit at which point they decay. The 
evolution of the oscillon as it approaches uj cr it is observed to be rather gentle, slowing down 
because the radiation of a quasi-breather decreases as we approach uj cr it from below. Fig. 
|3| reveals another possibility for oscillons, namely those starting from uj > uj cr it, in this case 
the quasi-breather picture would say that u approaches uj cr u from above with the oscillon 
following down one of the curves in Fig. ||. Using Gaussian initial conditions we were 
not able to find an oscillon with frequency larger than the critical value, however we are 
able to start with a quasi-breather profile satisfying uj > oj cr it- The results are somewhat 
disappointing for the quasi-breather picture, as the subsequent evolution did not show any 
oscillon phase. Instead, with the exception of D = 3, the field dissipated rapidly straight 
into radiation, uj = 1, as shown in Fig. |l^. 



5. Oscillon life-times and dimensionality 
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Figure 13: Left: The evolution of uj as predicted by the quasi-breather model using 4.7, Right: 
Lifetimes for an oscillon starting from a Gaussian with a = 3.3 for 2.5 < D < 7. 

In this section we shall collect together the information we have found about oscillons 
and quasi-breathers, and show how the quasi-breather picture predicts oscillon evolution 
should depend on dimension. 

We have already seen some specific examples of how the quasi-breather model predicts 
the evolution of uj using (|4.7|), in Fig. [l3| (left) we collect together the predictions for a 
range of dimensions 2^ < D > 5. As noted previously these tend to be off by a factor 
of a few, overestimating the lifetime, but we can clearly see that the oscillons in lower 
dimensions take much longer to reach the critical frequency. In higher dimensions we find 
that the predicted lifetime is simply not long enough for the system to enter into an oscillon 
phase. 

We may now compare the lifetime for some specific initial condition to that of the 
quasi-breather oscillon model and see how they compare in different dimensions 

As was shown in [?, ?], in D=3 one can tune a to produce oscillons with a large 
lifetime, but the level of fine-tuning necessary is extreme. It is then surprising to find that 
in D < 3, oscillons living for t ~ 10 6 are rather more generic and require little tuning of 
initial conditions. The converse is true in higher dimensions, with the amount of precision 
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required to find a long-lived oscillon increasing with D. This is the main reason that the 
results of [?, ?], based on substituting an oscillon profile into the action, needed to be 
confirmed with the more detailed analysis we have presented. Rather than trying to find 
the fine-tuned initial conditions which maximize the lifetime in each dimension, we shall 
take the same Gaussian initial conditions (<r = 3.3, A = 2) in the various dimensions and 
compare the lifetime of the subsequent oscillon with that of the quasi-breather model. The 
results are plotted in Fig. [l3| (right) and show that the semi-analytic predictions are well 
matched by the numerical evolution. 

In order to guide the eye we have overlaid a curve of the form 

tlifetime = (Arit - D)- 

with c = 1000, a = — 1.5 and -D C rit = 2.65, suggesting that the lifetime goes to infinity 
for some D > 0, and possibly for D > 2. While this is a possibility, the quasi-breather 
profile could conspire to have zero amplitude for the radiation modes, we believe that the 
extrapolation of the quasi-breather results to D = 2 require a finite lifetime. In particular, 
Fig. [?] does not suggest that the radiation goes to zero at D = 2. Note also that there is 
a kink in the curve of Fig. ^ (right) at D ~ 5, after which the lifetime suddenly drops 
(there are now only a handfull of periods), where a similar fit gives c = 45, a = —0.7 and 
-Dcrit = 5.02. Although the precise fitting parameters are not essential, the change in power 
suggests a different regime (no oscillons) has set in. 

6. Conclusion 

The interpretation of oscillons in terms of quasi-breathers gives a simple and consistent 
picture of why such long-lived objects should exist and may even be generated in a generic 
phase transition. 

In our setup, the quasi-breathers are the periodic solutions to the equations of motion 
with minimal energy. The presence of this radiation renders the quasi-breather unstable, 
although since the radiation is minimal, presumably it is the longest lived. Even without 
excessive fine-tuning, localised lumps (oscillons) can approach and follow the evolution of 
quasi-breathers, which hence must be considered (at least) weak attractors in the space of 
classical solutions. 

A quasi-breather exists for each choice of frequency. Given an initial choice of u>, 
the configuration evolves through a sequence of quasi-breathers, approaching a critical 
frequency, corresponding to the lowest energy quasi-breather. The time is takes to do 
so is well reproduced by considering the emitted radiation. As the critical frequency is 
reached, the energy can no longer be lowered, and the quasi-breather decays completely 
into radiation. 

We performed this analysis for a range of D, the number of spatial dimensions, in 
order to establish the stability/existence of oscillons in higher dimensions. Using the quasi- 
breather approach we have provided an understanding for the prescence of very long-lived 
solutions in low dimensions and short-lived solutions in higher dimensions. Although our 
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method could not cope with D < 2, it is tempting to conjecture from (|5.1| ) that the lifetime 
will diverge for some D > 0, although the extrapolation of Fig. |7] does not support this. 
Perhaps including more radiation modes and more hard work can settle this. While we 
have studied the dependence on D we have not studied the dependence on V(<p). One 
especially interesting case would be the sine-Gordon potential where we know there exist 
exactly periodic solutions in D = 1 [?] as well as D = 0, which is just a particle in a well. 

For D > 3, the lifetimes become smaller and smaller, and oscillons rush to reach the 
quasi-breather before decaying. For D > 5 we were unable to find quasi-breathers with 
our shooting method, and it indeed seems like a wholly new, oscillon-free, regime sets in. 
This is where the time a configuration takes to reach a quasi-breather is longer than the 
lifetime of the quasi-breather itself, so no oscillon can get established. 

In a 3+1-dimensional expanding Universe, assuming that the expansion does not sig- 
nificantly influence the field evolution, it is therefore natural to expect oscillons to be 
generated in phase transitions or during (p)reheating after inflation, with lifetimes of order 
10 4 in mass units, but probably not much longer. Scalar fields thermalise on much longer 
timescales than this [?, ?], but when including gauge degrees of freedom equilibraton is 
much faster (see for instance [?, ?]). Oscillons may potentially delay this process if they 
are produced in large numbers, and if the effects of the gauge field do not substantially 
change the picture. 

In order to make further statements about this point, we need to look for quasi- 
breathers in more realistic settings such as the Standard Model and its extensions. In the 
restriction to the 5'f7(2)-Higgs model [?], an approach similar to the one employed here 
should work, in spite of the complexity introduced by the additional degrees of freedom. 
Work is under way to address this. 
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